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Die folgenden Angaben sind den vom Anmelder eingereichten Unterlagen entnommen 

Priifungsantrag gem. § 44 PatG ist gestellt 

(54) Computertomograph mit reduzierter Dosisbelastung bzw. reduziertem Bildpunktrauschen 

(57) Die vorliegende Erfindung schafft einen Computerto- 
mograph (Ein- Oder Mehrzeiler mit oder ohne Spiralabta- 
stung), bei dem das Rauschniveau der interpolierten Pro- 
jektionen einen bestimmten Schwellwert nicht ubersteigt 
und der bei gleichbleibender Bildqualitat eine signifikante 
Reduzierung der Patientendosis eriaubt oder dement- 
sprechend bei gleichbleibender Patientendosis eine signi- 
fikant hohere Bildqualitat erzielt. 

Die zu rekonstruierenden planaren Projektionsdaten war- 
den nach folgender Forme! berechnet: 

die Filterung findet also entweder vor oder nach Umrech- 
nung auf Parallelgeometrie statt und minimiert somit den 
■ glattungsbedingten Verlust an Ortsauflosung. 
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Beschreibung 

Die modeme Computertomographie (CT) beinhaltet eine Vielzahl von MeBmodalitaten. Einerseits existiert die kon- 
ventionelle CT, also die Aufnahme einzelner Schichten. Diese wurde groBtenteils durch die Volumenaufnahmetechnik 
5 Spiral-CT ersetzt. Andererseits werden neuerdings CT-Gerate gefertigt, die statt nur einer Detektorzeile mehrere Zeilen 
(N Zeilen) besitzen. Mit diesen Mehrzeilencomputertomographen sind sowohl Schichtaufnahmen von nunmehr N 
Schichten gleichzeitig, als auch Spiralaufnahmen moglich. 

All diesen Modalitaten gemeinsam ist, daB eine Konvention ("z-Interpolation") X besteht, die gemessenen Daten p(p, 
a, n) (p ist der Winkel im Facher, a der Projektionswinkel also die Winkelstellung der Rohre bzw. Winkelstellung von 

10 Rohre und Detektorsystem und n zahlt die Nummer der Detektorzeile oder die Nummer der Schichtaufnahme) auf einen 
planaren transaxialen an der z-Position Zr giiltigen Datensatz px(P, OC, Zr) (Schichtdatensatz) umzurechnen. Dieser 
Schichtdatensatz wird dann mit einem Rekonstruktionsverfahren fur planare Daten (iiblicherweise gefilterte Riickprojek- 
tion oder Fourierrekonstruktion) rekonstruiert, um das CT-Bild zu erhalten. 

Bei konventionellen Aufnahmen (Schichtaufnahmen), die aus Einzelumlaufen an Positionen z', z", z'", . . . entstehen, 

15 wird iiblicherweise aber nicht notwendigerweise die Rekonstruktion nur an genau diesen Positionen zugelassen (d. h. zr 
e [z\ z", z"', . . . }). Bei Spiralaufnahmen ist die Rekonstruktionsposition Zr im allgemeinen frei und retrospektive wahl- 
bar. Die am weitesten verbreitete Konvention zur z-Interpolation bei Einzeilen-Spiral-CT ist der Algorithmus X = 
180°LI, eine lineare Interpolation zwischen in z-Richtung gemessenen Datenpunkten [Willi A. Kalender Wolfgang 
Seissler, Ernst Klotz und Peter Vock, "Spiral volumetric CT with single-breath-hold technique, continuous transport, and 

20 continuous scanner rotation". Radiology 176 (1), S. 181-183, Juli 1990], [Arkadiusz Polacin, Willi A. Kalender und Guy 
Marchal, "Evaluation of section sensitivity profiles and image noise in spiral CT", Radiology 185 (1), S. 29-35, Oktober 
1992]. 

Bei alien CT-Bildem gilt, daB die Bildqualitat in bezug auf Rauschen und Niedrigkontrasterkennbarkeit monoton mit 
der Patientendosis steigt. Das heiBt, daB der Anwender den Rohrenstrom und somit die Dosisbelastung des Patienten so- 

25 weit erhoht, bis die Bildqualitat seinem (subjektiven) Empfinden nach "gut" ist. Dosisbetrachtungen haben insbesondere 
in europaischen Landem einen hohen Stellenwert und somit werden viele Versuche unternommen, die Patientendosis bei 
moglichst gleichbleibender Bildqualitat zu reduzieren. 

Eine einfache Moglichkeit zur Reduzierung des Bildpunktrauschens bietet sich durch die Wahl eines glattenden Re- 
konstruktionsfilters. Die Rekonstruktionsfilter konnen an jedem kommerziell erwerblichen CT-Gerat in gewissen Gren- 

30 zen frei gewahlt werden und somit kann jeder Anwender das Rauschniveau im Bild verringern, ohne die Patientendosis 
zu erhohen. Da mit dieser Methode der gesamte Datensatz geglattet wird, geht dies zwangslaufig mit einer Verschlech- 
terung der Ortsauflosung einher. Das Problem, die Patientendosis bei moglichst gleichbleibender Bildqualitat zu reduzie- 
ren, bzw. die Bildqualitat bei gleicher Patientendosis zu erhohen ist damit nicht gelost, denn es muB ein KompromiB zwi- 
schen Ortsauflosung und Bildpunktrauschen gemacht werden. 

35 Zur Losung des Problems finden sich in der Literatur Ansatze zur adaptiven Filterung der MeBdaten, d. h. der Daten- 
satz wird nicht global, sondem nur lokal geglattet [Jiang Hsieh. "Adaptive trimmed mean for computed tomography 
image reconstruction", Proc. of SPIE, 2299, S. 316-324, 1994], [Jiang Hsieh, "Generalized adaptive median filter and 
their application in computed tomography", Proc. of SPIE, 2298, S. 662-672, 1994], [Jiang Hsieh, "Adaptive filtering 
approach to the streaking artifact reduction due to x-ray photon starvation". Radiology 205 (P), S. 391, 1997], [Berkman 

40 Sahiner and Andrew E. Yagle, "Reconstruction from projections under timefrequency constraints", IEEE Transactions 
on Medical Imaging, 14 (2), S. 193-204, Juni 1995]. Ublicherweise werden zur adaptiven Filterung benachbarte Detek- 
torelemente verwendet die Filterung findet also ausschlieBlich in p-Richtung statt. 

Ansatze, die Filterung rohdatenbasiert adaptiv (angepaBt an jeden einzelnen gemessenen Projektionswert) in alien drei 
Dimensionen P-, a- und z-Richtung) durchzufiihren sind bisher nicht bekannt. 

45 Der in Patentanspruch 1 angegebenen Erfindung liegt die Aufgabe zugrunde, fiir einen Computertomographen (Ein- 
oder Mehrzeilersystem mit oder ohne Spiralabtastung) das Bildrauschen und die storenden Rauschstrukturen, die durch 
Korrelation benachbarter Bildpunkte entstehen, weitestgehend zu reduzieren. Damit wird eine Verbesserung der Bild- 
qualitat bei gleichbleibender Patientendosis erreicht bzw. eine Dosisreduzierung bei gleichbleibender Bildqualitat mog- 
lich. 

50 Diese Aufgabe ist erfindungsgemaB gelost durch die Merkmale des Patentanspruchs. 
Die Erfindung ist nachfolgend erlautert. 

Es bezeichne px(^, Zr) die planaren Projektionsdaten (z. B. Schwachungswerte) eines CT-Scans, die an der z-Posi- 
tion Zr giiltig sind. Dabei ist ^ der Kanalindex und entspricht in Fachergeometrie dem Winkel im Facher und in Parallel- 
geometric dem Abstand des Strahls zum Rotationszentrum. ^ ist der Projektionsindex und entspricht sowohl in Facher- 

55 als auch in Parallelstrahlgeometrie der Winkelstellung der Rohre bzw. der Winkelstellung von Rohre und Detektor. X sei 
der Name der Konvention zur z-Interpolation, mit der diese planaren Daten aus der Messung erzeugt wurden. Die Re- 
konstruktionsposition zr sei retrospektiv auf einem der Messung entsprechenden Intervall [Zmin, Zmax] fr^i wahlbar. Dazu 
drei Beispiele. Beispiel 1: Im Falle der Einzeiler-Spiral-CT ist X = 180°LI ein denkbarer z-Interpolations algorithmus, der 
hier zur Anwendung kommen kann. Beispiel 2: Bei Mehrzeiler-Spiraldaten kann, analog zu Beispiel 1, eine Standard z- 

60 Interpolation zugrunde liegen. Beispiel 3: Bei konventioneller CT (Schichtaufnahmen) kann sowohl bei Ein- als auch bei 
Mehrzeilem als Konvention X eine Nachster-Nachbar-Interpolation oder eine lineare Interpolation zweier benachbarter, 
gemessener Schichten verwendet werden. 

Bei dem hier beschriebenen Rekonstruktionsverfahren fiir die Computertomographie werden die planaren Projekti- 
onsdaten px(^, Zr) vor der Rekonstruktion adaptiv nach folgender Formel gefiltert: 

65 

PAF(t ^, Z) = Jd^'dd'dzTA^(^-4')fA^(^-^')fAz(z-z')px(^', ^\ z'). 

Dabei bezeichnen 
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- px(5, z) die vor Durchfiihrung der adaptiven Filterung (AF) zur Verfiigung stehenden Projektionsdaten in Par- 
allel- Oder Facherstrahlgeometrie. Diese errechnen sich aus den MeBdaten durch die Konvention X ^ AF fiir belie- 
bige z. 

- Paf(^, z) die nach Durchfiihrung der adaptiven Filterung zur Verfiigung stehenden Projektionsdaten in Paral- 
lel- Oder Facherstrahlgeometrie. Diese Daten konnen dann der Rekonstruktion (z. B. gefilterte Riickprojektion) zu- 
gefiihrt werden. 

- A^, A'd und Az die Filterweiten in den drei Koordinatenrichtungen. Die Filterweiten sind Funktionen des aktuell 
adaptiv zu filternden Projektionswerts px(^, z): A^ = A^(px(^, z)), A'd = A'd(px(^, i^, z)) und Az = Az(px(5, 
z)). 

- fA^( • )j fAd( • ) und fAz( • ) die Filterfunktionen (achsensymmetrisch mit Werten > 0 und Gesamtflache 1) fiir die 
Glattung in den jeweiligen Koordinaten. Die Filterweiten A^, A'd und Az stehen jeweils fiir die Halbwertsbreiten 
Oder ein anderes charakteristisches BreitemaB der Filterfunktionen. Wird eine oder mehrere der Weiten gleich Null, 
so reduziert sich die Filterfunktion zu einer Diracschen Deltafunktion und es findet in den entsprechenden Koordi- 
naten keine Filterung statt. 

In einer moglichen Ausfiihrung konnte die Wahl der Filterweiten mit der Formel 



(A^ + #)(A5 + i«)(Az + 2) = c,^^ 



gPx(?,».^) falls g/'^«.f-:^)> 7-2 
sonst 
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erfolgen. Dabei bezeichnen 

- ^, 'd und z den Abtastab stand in den drei Koordinatenrichtungen: ^ ist die Breite eines Detektorelements, ^ ent- 
spricht der Detektorintegrationszeit und z der Halbwertsbreite des Schichtempfindlichkeitsprofils. 

- e^^^^' ^' die relative Varianz des Projektionswerts px(^, z) als MaB fiir das Rauschniveau des Punkts z). 

- den Schwellwert ab dem die Filterweiten A^, A'd und Az einen von Null verschiedenen (positiven) Wert an- 
nehmen. 

- c die Normierungskonstante, die so gewahlt wird, daB bei Unterschreitung des Schwellwerts T die Filterweiten 
A^, Ai& und Az den Wert Null annehmen. 

Folglich gilt cT^ = _ 

Die Bedingung (A^ + ^) : (Ai^ + : (Az + z) = ^ : : z stellt sicher, daB das Verhaltnis der effektiven Filterweiten (dies 
sind die Filterweiten plus der intrinsischen Mittelung durch das CT-Gerat) dem geratebedingten und somit vom Herstel- 
ler optimierten Verhaltnis aus Detektorbreite, Integrationszeit und z-Ausdehnung entspricht. 

Die Filterfunktion selbst kann z. B. eine der folgenden Darstellungen haben: 

1 fl falls 21x1 < Ax 
• Rechteck: fM = — \ 

Ax 0 sonst 
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Dreieck: 



1 fl - \x\/ Ax falls Ixl < Ax 
Ax 0 sonst 
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Gauss: 



1 



2^ Ax 
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Diese Darstellungen haben den Vorteil, daB die zugehorigen Integrationen analytisch ausgefiihrt werden konnen und 
die adaptive Filterung sich auf eine gewichtete Sunmie uber die gemessenen Daten reduziert. Damit lassen sich die ad- 
aptiven Filter effizient implementieren. 

Patentanspriiche 

1. Computertomograph (Ein- oder Mehrzeiler mit oder ohne Spiralabtastung), bei dem durch geeignete Filterung 
sichergestellt wird, daB das Rauschniveau der Projektionen einen bestimmten Schwellwert nicht iibersteigt. 

2. Computertomograph nach Anspruch 1, bei dem die Interpolation der Projektionen adaptiv erfolgt (adaptive Fil- 
terung, adaptive z-Interpolation). 

3. Computertomograph nach Anspruch 2, bei dem die Projektionsdaten in Parallel- oder Facherstrahlgeometrie so- 
wohl in Kanalrichtung (^-Koordinate), in Projektionsrichtung ('d-Koordinate) als auch in axialer Richtung (z-Koor- 
dinate) adaptiv gefiltert werden. 

4. Computertomograph nach Anspruch 3, bei dem die Projektionen adaptiv nach folgender Formel gefiltert werden 
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(3D- adaptive Filterung): 

PAF(i ^, z) = Jd4'dd'dzTA5(^-^')4^(i!>-d')fAz(z-z')Px(^', z'). 

5 5. Computertomograph nach Anspruch 4, der bei nahezu gleichbleibender Bildqualitat eine signifikante Reduzie- 

rung der Patientendosis erlaubt, oder dementsprechend bei gleichbleibender Patientendosis eine signifikant hohere 
Bildqualitat erzielt. 
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